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GOMPRjESSLBLE POTiSNTIAL ELOW ■ WITH-.-.CIRCULATION 
■ ■ .ABQ.UT a;. CIRCULAR CYLINDER " , / ., 
By Max A, Heaslet . 

; : < . .SUMMAHY 

The potential function "'for flow, witii' cir culat ion^ 
of a compr es s ilDl e fluid . alD out a circular cylinder is ob- 
tained in series form including t erms- 'of . . t he order of M"^ 
where M is the. Mach numher of the free .streaiRo ' The re-^ 
suiting equations are used to ott ai n pr es sur e coefficient 
as a function of Mach number at a point o,n the surface of 
the cylinder for different values -of circulation^. The 
coefficients deriv.ed are ^compared with- the &lauer t-Prandt 1 
and Kar man-Ts i en approximations v/hich are functions of the 
pressure' coefficient s of an in'compr ess ihl'e fluid. For the 
cases considered;, the values of the pressure coefficients 
computed from the theory were 'found to.. lie somewhere be- 
tween the two approximations J the. first underestimating 
and the second overestimating the theoretical result, 

INTRODUCTION 

In the tw6-dime nsional irrotatipnal flow of a com-- 
pressible fluids where the expansion is assumed to be 
adiabatic'5 the velocity potent ial is known to satisfy a 
nonlinear partial differential equation of the second 
order. For subs^onic velocities, at least three methods 
are known for the approximate solution of. this equationo 
They are usually denoted as the method of small perturba- 
tions p the Rayleigh^Janzen method^ and the hodograph 
method. 

The method, of small perturbations (references 1 and 
2) assumes that velocity changes which are brought about 
by the airfoil in- the uniform parallel air stream are 
Bmall in comparison with .the velocity of .the undisturbed 
stream* Under this assumption it is possible to intro«. 
duc3 new variables which reduce the differential equation 
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to a Laplace equat i on, -and, as- a ^co.ns e.q.ue the problem 

oecomes one concerning flow in an incorapressil^le fluid, 
provided the TDody is assumed distort ed to correspond to 
the change; of variahles.- ; The • assumed distortion . cons.ists 
in expansion of the dimensions of the airfoil perpendic- 
ular to the direct ion of the fre-e .stream in the ratio 

1/ \-l^K , v^here M is the Mach number .of the undisturbed 
stream* 

The Eayleigh- Jansen method (references 3 ^-^^ ^) 
sumes that, the general expression for velocity potential 
may "be written as a series in rising powers of M and 
with .variable .cp efficients . .These • coef f icif nt s .can. he 
shown :t 0 . sat isfy ^c er tain ,Poissoh differential equations 
andy if the equations; are int egr abl e the - so lut ion,""h e- ^ 
comes a-, matt er :0f -.det.erm.ining these co ^f f i ci ent s « Sue-, 
cessive steps, .however, /become increasingly laborious and. '.', 
the convergence . of the .-series . may be sloWg even at. rela- 
tively small -Mach-numbers J if -the - shape of the bodyis 
such that the.- speed of sound, is approached locallyo • . Solu- 
t i ons , .using- this method of attack^ have b een carri.ed out . . 
by Oo v^Kapian • ( r ef er ences - 5- and 6) , S, Gc ' Hoqker' (reference 
7 ) 3 .1 o ■ I m a i (references). K p Tama d a^_an d„,Xo-:^c a-H'0--^t . ' 

encQ, 9) , and X. ..Poggi-',('r"ef'^^ Poggi introduced 

certain r ef inem.ent s , > and some of the ..p^ ec eding ref er en'c es 
employ this processr It ig tantamount .to using .the so- - 
called Neumann function in solving given Pois son " equa- 
tions and will be discussed latere 

The hodograph method is ascribed by writers on that 
subje.ct- to Pe Molenbrock and- Ac, Ts chaplygin^. Instead of 
express ing the , velocity pot ent ial as a function. of coor- 
dinates in the . Car t e^ian or polar plane. th6 magnitude., 
of velocity V . • and its inclination 9 r t 0 . an , as sumed 
axis ■ ar? ' chosen as > independent variables The resulting- 
differential equation, is linear and can b e f ur ther . s impli** 
fied by replacing the- pr es sure-volume relationship f.or 
adiabatic expansion by • the • equat ion of a line tangent at . 
a point corresponding to the state of the fluid in the 
ambient streamc This artifice was suggested by To von 
Karmdn (ref ere-iices . 2 and 1.1.) .. and us ed mos t., successfully 
by Ho. To Tsien ( refer ence.;12) K* Tamada (reference I3) . 
also has applied Ts i en^s • more general result s on., ^llipt ic 
cylinders to compressible flow past a circular cylinder o. 

• One. nct.eworthy. result of the. hodograph method has 
been the K4rman-*Ts i en expression for pressure coefficient 



3 



P in .terms of Mac h . number M- and- Pj.(-.q» ' tho pressure 
coefficient for -U-cy:' S-his "expression m-ay 'be v/ritten 



•F Pm=.o ^-^ — — • (1) 



s 2 



It alwaj^s- gives for. nogati -re:, p.ressure; coefficient s > a 
result great'Or -.in- absolute value titan the (xlauert^Pr andt 1 
formula^ whlcTa is: based on the method of small perturba- 
tions,- 



•• ■■■■ ■ -P = Pm=o 7-::i::i:r- • (2) 

and is currently accopted as the more accurate of the tv;o« 

From equations (l) and (2) it is possi"ble to com- 
pute the critical Mach number Mq, the value of M . at 
which the local speed ' of sound is attained, in terms of 
P|j_Q .• The relations involving Uq and Pm=o» corre- 
sponding respectively to formulas (l) and (2), are 



A. f(.A..y^M^- -1] = .^., 

YM^ i\Y+l Y+1 / J / Ut P,, 

<V ■^-^''c + . 2 

l+n/ 1-M^ 
P 

and 



(3) 



%=o = Ml - (1 - (4) 



The difficulties inherent in the last two procedures 
are quite as distinctive as their respective approaches 
to the prohlom. As stated "before, the Rayl eigh- Janz cn 



iiiethod employs classical mathemat ics ^ ■ the required terms 
teing solutions o.i Poisson equations with given "boundary 
conditionsj tut the work involved is arduous* In the 
hodograph met hod '"the •.principal difficulty is to deter-- 
mine proper "boundary conditions in the Y^Q plane« In 
available calculations the solution is given with a 
slight distortion in the given houndaryo It is possible 
to correct ..this-. distortion, in some - cases so that the 
finsLl r 6sult s ^.a:re not ' too --s eriously .-affected* *■ When the 
flow "a]i:*ound '.t3ae' t cir culat ion^, howeverp the 

change in the boundary is more serious^ for nonperiodi.c 
terms appear and the boundary is no longer a closed curve* 
At the present time, no way has been found to circumvent 
this trouhle^ Added cir-cul-at ion does not involve any 
essential variations in the. Rayleigh-Jans en method, how- 
ever, and in this report the velocity potential for such 
compr es silDle flow ahout a circular cylinder has "been de-. 
riy.eda ,-Sin.ce no theoretical study has teen presented, as 
far as is known, to determine the error in the Karman- 
Tsien pressure: coeff icieht,' the' results ohtained in this 
report^ fur nrsh ra means -of a.pproachirig this protleni. The 
results of... such;; caTcuiat i-ons., f'O'r various values of ...cir- 
culat ipn.p=: a-r e |jh,er ef o;r e . in;Clud:ed.V ■ ■ ' ' ■ 



AIALISIS 

Consider a gas obeying the adi^tatic law. :;and -.flowing 
irrot at ionally in tv/o dimensions. Its equation. 'of motion 
may he written in- 'polar .coordinat es in the form 

Ill ui^l ^•■i\y$ = L.'^'(B il (5) 

where 

$ velocity potential 

3_ ■ 3^ _^ -i- 80' -'i ■ 33$ ', ;', . ' .' . 

8r^ : T - Bf '-i^ bp 
V ratio of specific heats of gas 



velocity of sound in undisturted. flow 
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U velocity .of free stream 

M = — Mao.h num"ber of free stream 
°o 

.^2 ^ fb^ Y ^ JL^ Y local velocity squared 

W J Ue J 

By introduction of the varia'bles cp and v, so that 

$ ' V 

cp = '-> and V = -> 
U U 

equation (5) may "be v;ritten in the form 
'Y-1 



1 - 

2 



M= (,=-1)1 v=9 = i >«= kT^'hTT^) 



v/her e 

.2 



.2 ^ /S'P 



/in 1- /.^ie\ (7) 



following the method of P.ayleigh and Janzen, assume that 

cp may "be developed in a series of ascending powers of 
M^j so that 

cp = + M^cpi + M^cpg + • o . (S) 

After substitution of equation (S) in equation (7), ele- 
mentary calculations shov/ that 

«'^- 0 1 « 

where 



° Ur ; r^ / 



(10a) 
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= 2 i 



1 



2 ^ "^5 



(10b) 



(10c) 



• « 9 O 0 



In a similar manner, equations (S) and (9) may b e - s'libst i- 
tuted in equation (6) and on eqviatin^ coefficients of the 
same, powers ..of .H., ..the follpv/ing .r.ela.t ions for.. 9g , 
^g, . . o result: 



= 0 



9, - 



3 



1 



2 \Sr br he bQ J 



(11a) 
(lib) 



= i (Y ^ 1) (v^ - a) ^^cp ■ 
2 0 * 

+ 1 (^^..^^ J. ^ ^) 
.. :2 \9r : . r^ ,^8 ^6 /• 



(11c) 



9 • C 



If eciuations (lla), (lib), and (11c) can *be solved 
successively for 9o,9i9^33 • • » the values may be sub- 
stituted in equation (2) to get the potential function for 
the flow of a compressible fluid. A s t ep~.by-st ep proce- 
du.re is therefore established v/herebj" any desired degree 
of approximation to cp ' may ^be . obtained, ..provided the 
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value of M is v;ithin the region of convergence of the 
resulting series^ ^ Equation (lla) is the differential 
equation satisfied loy the potential function in the case 
of incompr essitil it y» Once this potential function is 
knov/n, it is used to evaluate the rig:ht-*hand memlDer of 
the second equation, the solution of which furnishes the 
second term in the development of cpc The method of oh« 
taining further terms ^follov/s the saine general procedure. 

Consider now the case of a right circular cylinder 
of infinite length in a compr es s ihl e fluid, the axis of 
the cylinder hei'ng'at right angles to the ^direction of 
steady flow« In determining the velocity distr il^ut ion 
atout the cylinder, the prbhlem may he treated two- 
dimens ionally with a circle as the houncijary curve and the 
equations established in the Hayleigh- Janz en method may 
he applied directly, in the following manner* The radius 
of the circle i^ arhitrari.ly assumed' equal, .to 1, and a 
polar coordinat-e system is chosen with origin at the cen- 
ter of the circle and polar axis extending down streamo 
The flow ahout the circle is assumed to he that resulting 
from the comhination of uniform stream velocity and cir- 
culati-on ahout the .cylinder . Under these, conditions, the 
classical expression for cp^ is \\rel*r knowh^ It may he 
writ t en 

/ 1 , * r 

cp = (r 4- J cos ' 0 « — 0; 
.0 V r/ 2ttU 

where F is the circulation around the circle, imeastired 
positive in a clockwise direction. For ease of computa- 
1;ion it is convenient to set 



ttU 

and, as a consequence, 



/ 1 \ K 



r ,/ .-.2 
The "boundary conditions, in general, are 



g = 0 for r = 1 (13a) 



and 



^ = cos 0 ' for r = CO (13"^) 



From equations (l2) and (lOa) 
^ = ^1 - ^ cos PG + K^i + Ljsin 9 + (lU) 
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This ,result, together with eqviation (llh), gives 



= ( :3 :?] cos 6 + ^ cos 3 8 . 



V r= rV 

(i- 1 ' 2 1 \ K ^ 

« + ^—1^) u — COS 9 (15) 

2r^ 2r^ 2rV 2r^ 

The more elementary methods of integration lead to 
certain difficulties v;hen an attempt is made t o solve f or _ 
cp^ 5 in equa t ijCLn_(_l-5-)-» T-h e s-e— d-i-f-f^rcuirt i"e"s FeiHl t from 

nonperiodic terms in the particular integral and resultant 
trouhle in determining such constants of integration that 
tho necessary periodicity^ in terms of 9^ is maintained 
in the final expression for the potential function. This 
difficulty may he ohviated, hov/ever^ hy established meth-- 
ods. (See appendiXc) It follows that the solution of 

2^ sinmQ /-/-s 
V Q = (16) 



satisfying the boundary conditions 



1 V / r=co 



IS 



sin m 8 



jii(in-s + 2) (m+s-2) 

and 



ilzil « ^iL. I (when m^2h) (l7a) 



Q = ZllRJLl fl. + I log rl (when in+2 = s) (1?^) 

mr^ lj.m- 2 J. 

The veracity of these soltitions, together v;ith analogous 
ones existing when -sin m 8 is replaced "by cos m 9, may 
"be checked easil;?* TDy suhstitut ion in eq.uation (l6). 

Since eciiiation (I5) is a linear differential eqiia- 
tion, its solution is determined "by considering!; each term 
of the r.i^:ht-hand memlDer and summing the individual inte- 
f^rals ohtained hy means of equations (l7a) and (17'b), 
She final result is 



= cos 9 

1 



ii- - Jz^ + — "^1 + cos 3 9 
I2r 2r I2r / 



\hr 12r7 

'-') 



+ Ksin29(^ + -i- + --^ + 

VS 6r^ 2^4t'^ 2r^ 



+ cos e (.— + ) (12) 



\hr hr J 

In the evaluation of cp the calculation follov/s 

2 

the sam.e pattern of development. Prom equation (lOh)^ 
together with equations- (12) and (IS) 

+ cos U 0 ( — I + K sin 9 j — + — — ^ + Ji~ + L 

\vy Ur 6r=^ I2r^ 3r' 



log r \ 



10 



+ Z sin 



3 e f'_ A - -1. -i- + - 2 log A 

\ hr 6r=^ 2r^ I2r' ; / 



+ I + + 

\Ur^ Ur^ 2r^ / 

+ cos 2 e f-nl - ^ - -J- - 
+ K^- sin e f + + ) (19) 



1 1 log r" 





I2r* I2r^ 






■ 1 


log r"^ 







This result, together v/ith equations (l2), (l^), and (IS), 
substituted in equation (lie),' gives 



- 7~3 l^ ' 

+ cos 3 9 — - I + cos 5 9 — 



+ Zsin2 eY— L^-J_ + - 3 



2r* Ur® 2r^ ^^r^" 



+ Z sin U e + 



^2 . /' 3 9 3 
+ K*" cos 6 - - + — 

V2r.=* Sr^ gr'' Sr 



9 



+■ cos 3 S — r- + •+ - 

\gr=* gr^ gr' gr^- 
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+ sin 2ef-^^^ ^ + + cos 9 ( 

VlGr"* gr^ l6r«/ U6r7j 

r 

! /-32 39 15 11 \ «/ 19 3 5 1 

+ <^ cos Si -^+-^ + ]+ cos 3 9/ — ^+ 

l_ V3r^ ar*^ 3r^V -.• - \Sr^ 3r^ 2r^ 



+ cos 5 e(.^ - -i^) ■ 



. ^/ -l 11 .25 1^+^ 3 6 log r 4 log 

+ K sin 2 9| : — ^ ;+ =-r - -*^+— =^ - ^ — + 

' ykr" 3r* I2r® 3r^ 2r^o 



, / 1 11 ? 5 1 3 log r' 
+ K sin 4 9 + -t -» — ^+ • — + . ■>■■ ■■ . 



/~3 2 13 log r log 

+ K cos 9 - — — — - — + — 

\Sr3 3rS Z^r"^ 2r^ , 



,3 , ^3 ^ 2 ^ 3 1 ^ log r ^ 3 log r log r 



+ Y." cos 3 e| 



\Sr' hr"^ kr^ 2r 



r r 



+ sin 2 9/^^ - -i. - + k4 cos ei^.::! 

V2r4 3j.e 2r* / U6rV J 



(20) 



To integrate, formulas (17^'^') and (1713) are again 
resorted to. The method of integration given in the. 
appendix also provides integrals cor r esponding , t o the 
nev/ type of terms appearing in the right-hand member of 
equation (20)« Thus, the solution of 
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. i£i.£4i^_^Lj (21) 

satisfying "boundary conditions Ur*- " ^» Is"^ 



IS 



sin m e r 1 ^ 1 1 sin m 9 | 2- m log r 

^" " "2m"~ |(in4-s-2)3 Tm-s + 2)2j 2inrS-2 1 (m+s-2) (in-s + 2) 

+ 3: - ^ I (22a) 

(m+s-2)2 + { 

when m s-2» When m = s»-2s 

/ 



-sinm9/l- o .1, 1\ \ 

... log^ r + — log r + —t] (2213) 

2m \2 2m J 



Proceeding di re^ ljv^_w^i t h t h o i nt_ef:r-a-ta^o-n --r-e-a-u-l-t-s — i-n— t-h'S" 
following expression 

CP, = (V-1) [cos 9 - * ^ . ^ . 1 ' 



3 <j^COS ^ gQj! i2j,5- ]^gy7 goi.% 

/ -61 3 1 \ 

+ cos 3 9 — *. — + — - 1 

■ V2>+Cr= l6r^ UOrV 

+ cos 5 S - -~^\+ K sin 2 9 " T^i 

■ ■■ _ ..,3.,„ + ^ ' + log r \ • 

6Ur® gOr^ Sr2 i 



/ -1 1 13 log r' 

K sin U S I + 



\l6r^ 320r* gOr® l6] 



4 



p / ^9 3 3 1 log r 

cos 8 — + — + — — 

Vl2gr 64r=^ G^+r^ I2gr^ l6r 



K= cos 3 ef-' " ^ ^ ' 3 



sin 2 e 



\6Ur GUOr^ I2gr^ 320r'' iGr^ 



^1 ^ _1_ ^ 1 + 3 log r \ 
V512P^ 32r* 5l2r® GUr^ / 



K"* cos 6 



+ 3 1 



I2gr I28r' 



'J 



< COS 9 I *- —is + — — = « + .... .u 

I VsOr 3r= l6r^ l6r' 240p 



) 



COS 3 9 



;:i:5 „ X. ^ _3 . 

,^6r Ugr=^ l6r^ 2^r'' , 



/I 1 1 

cos 5 9 / ---- + 



\^l6r I6r3 aOr^/ 

( 1 2267 23 19 1 

K sin 2 9—5 + ■ ^ - . ■ . ■ v x - — ~ + r- — ?, 

V16 2gg0r^ 144r^ 192r^ 40r^ 

+ 11 log r ^ log r ^ log r \ 
I2r3 "* 2r4 gr/ j 



Ik 



k / ~* 7 7 - 1 1 log r log r\ 




. cos 9 (2^ . JlI. .. Jl^ + 12^^ - + ilO 



\2ggr 96r3 7.2r^ l6r Sr=* 



/ -1 , 911 1 1 log r 5 log 



+ cos 3 e(-"r- - — r* + 



Vl6r 5760r=5 IgSr^ iGOr^. . l6r 12r=* 

) 



log r log^ r' 



r,3 . o / ' ' -■ . 5 log r log^ r^ 

+ sin 2 0 [ —"--^ — + ~ — + r- 1 

32r^ l6r^ / 



+ K'^ cos 6 U-i- ~ — i-g-^i (23) 
Vl2gr i2Sr"'y 



Applications of Theory 

Vifith the expreasions for cp^, (p^s and cp^, the two- 
term appi'oxima'tiwn for velocity potential is 

$ = U /' CP + <p. + cp M^^ izk) 

1 3 / 



and from this function the values of velocity at any point 
in the plane may be computedo Of particular interest is 
the evaluation of 

1 0$ 
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for this gives velocity normal to the radius vector of 
the point in question and thus, when r = 1, is equal to 
the velocity at the surface of the cylindert 

Neglecting all terms containing powers of l/ r 
greater than the first, Glauert (r ef er ence l) has given 
this well-known result 

1 ^^^^ . u /gin 6 + i- -^£ll^l— ] 

and under the same restrictions equation (sU) gives 

1 L*^ = U ('sin 9 + ^ - Kjf^os^ - KM^cqP + EtB2I 
\ T hQ J \ 2r ' 8r l6r 



These results are identical to the order of M <, 
Velocity at the surface of the cylinder is 



V (1„ G) = U 2 sin e + I + M' 



2 sin 9 sin 3 9 
5- I 



2 K cos 29 K^sin 9 

3— * -IT- 



+ M4 



(■Y-l) I ^ sin 9 
\120 



11 . , « 1 ^ A 127 „ 

; — sin 39 + - sin 5 9 - --r- K cos 2 9 

Ho g 2^0 



+ |J K cos H 9 + 1^ sin 9 - sin 3 9 



SO 



320 



■"■•^ra K^' cos 2 9 + ^ K* sin 9^ + sin 9 



12s 
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25 3 887 

i sin 3 e + - sin 5 9 - K cos 2 6 

2l+ g 720 



+ 55 K cos U e + — sin e - — sin 3 9 

72 lUU 960 



61' 
288 



cos 2 6+ 



1 

Pi 



s m 



(25) 



Results derived from this equation will te presented 
in terms of pressure coefficient, which is defined as 



(26) 



where the zero sutscripts refer to free-^stream conditions. 
By means of Bernoulli's equation it follows that 



M=o 



(27) 



where Pmsq cL^notes the pressure coefficient for incom- 

pressilile fluidso If P denotes pressure coefficient 
for a compr essible fluid obeying the adiabatic law^ then 

V 



P = 



r 



Y M' 



LL 



v-i , 

1 + M" 

2 



1 - 



1 



- 1 



•J 



(28) 



where M is the Mach numher of the free stream and 7 
is the ratio of specific heats (l.Ho for air)* 

As an approximation for P the Glauert-*Pr andt 1 re- 
sult (reference 2) is given "by equation (2) and the 
Karman-Tsien result (reference 2) is given equation- 
(!)• 
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The velocity at the topmost part of the cylinder may 
te found by setting 0 = 90^ in equation (25) and the 
resultant expression is a function of K and M. "In . 
figure 1, pressure coefficient at this point is plotted 
against M for K = 0, As a test for rapidity of con- 
vergence the expressions for velocity, using only M 
(one-term approximation) as well as and M' (two- 

term approximation), are used. It is to be noted that 
the curves diverge greatly near the critical Mach number, 
but that, for smaller values of M, the curves derived 
from equation (25) are together and definitely lie^between 
the results derived from, the Glauert-Pr andt 1 and Karman- 
Tsien relations. Figures 2 and 3 show the saDie equations 
applied for K = l/U and l/2v respectively. It. thus 
appears from these calculations that the true value of P' 
lies somewhere between the approximations applied, .On 
the other hand^ experimental data, as determined frorji air- 
foils, have shown a much better agreement with "the Kar'man- 
Tsien equation than have the theoretical results obtained 
here for the cylinder. 

In figures ^ and 5 the same point on the cylinder is 
under consideration, but circulation is made negative by 
setting K equal to --l/^ and -I/2 in the two cases. As 
the pressure coefficient gets smaller in absolute value, - 
the theoretical data agree more nearly with equation (l). 
For all the calculations, the one-term and two-term ap- 
proximations diverge \^idely as the Mach number increase's-. 
This is to be expected for, as has been pointed out by 
Messrsc, tte I, Taylor and C. Sharman in reference 14. 
the convergence of the series fails when -M reaches it.s 
critical value^, For near-critical velocities, several 
more terms would be required to furnish an accurate eval- 
uation of the true potential-flow pressure coefficient. 

Figure 6 shows the value of pressure coefficient at 
all points on the surface of the cylinder. These results 
were derived from equation (25) with K set equal to l/^ 
and at a Mach number of l/^*- Crosses on the graph are at 
positions obtained from equa'tion (l), and the circles 
were determined by equation (2)» The di sagreement " at the 
extreme pressure coefficients is again in evidence. 
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APPENDIX 



For the integration of the differential equation 



1 + i ^1 = f (r,9) 



Sr r br r 39 



(29) 



with "boundary conditions 




and 




= 0 



(30) 



it is' assumed, as in Poggi^s method, that the unit circle^ 
with center at the origin is surrounded externally "by a 
continuous dist ril:ution of sources such that the source 
strength of an element RdRdu) is 



Equation (29) may then te interpreted as the .equation for 
incompressible flow in such a region. 

The velocity potential of an incompr essi^ble fluid at 
point (rp9) due to a unit source at (H^u)) ma.y l^e calcu- 
lated "by the method of imageso If this potential is de^ 
noted "by cp, then 



f (E, 0))HdE"da) 



(31) 




i log r^ + - 2Rr cos (G - 0)) 



cp = 



2TT 1 2 




2 



r 



cos 




The required potential cpp satisfying equation (29), is 
therefore 
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cp(r,e) = 




[ 



log + - 2Rr cos (6 - u) ) 



] 



+ log 



~ + ~ 2? cos (e - w) 
r 



] 



+ 2 log - I f(R,«)) HdHdu) 
E 



(33) 



where the integration extends over the region of the 
plane lying external to the unit circle. 

In the eq^uations under consideration in this report 
the function f(E,a)) is restricted to one of the forms 



sin mo) 

M ■» ■ I I ■ 



COS mu) 



R 



s 



log R 



s m muj 



R^ 



cos mu) 
log R — ^ 



s ^ 1 

As an example of the integration process, take the first 
case listedo Then, set r^ = l/r, which results in 

. 2.7T 



cp(r,e) =.A P r log (v^ + R^ - 2Rr cos (8 - J)\ 

1 0 



sm mu) 

RdRdu) 



CO 3Tr 
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27T J J 
1 0 



log H — 5* Edadu) 



Ii + Is + I3 



Integrating I3 first v/ith respect to ' o) shov/s 
immGdiately that its valtie is zcro« ' : 



form 



For pnriposos soon evident is written in the 



1 0 



+ leg 



1 + / '^ly cos ( e - uj) 



s m m<jo 

'as 



adHduD (35) 



Since the log a - term vanishes, after integration with 
respect to the. expression for may "be simplified 

further hy the sulDst itution 

log 



■ • (0 



2~ cos (e - w) 

a 



n=Ti 

s ince r ' < a, and 



COS n (9— 0) ). 



(36) 



Ob 




TT 



I rz ,~ I I % ~ — ^ — — sin mtu (cos n 9 cos nu)- - 

^ .2TT J J n UaJ a""i 

1 0 • n=:i - 
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00 



+ sin nG sin no)) dRdu;= - ' ^ 

2tt 



m p ini[^m+ s - 1 



-TT sin.m 9 dP. 



1 sin in9 

2 mr^ m4-s-2 



(37) 



To integrate I^, the region exterior .to the unit 
circle is hroken into tv;o parts: The first part is a 
circular ring external to the unit circle and extends to 
the fixed point r; the second, part is the remaining por- 
tion of the plane and extends infinitely. Then 



I, = 



I J \ log r^^ log 

^10 L 



1 + 



2~ cos ' ( 9 U)) 
r 



sm mcu 



RdRdu) 



r 0 ^ 



log R"" + log 



1 + 



2(-l cos (9 - w) 



sm mu) 
... ■■■ k 



J 



RdRdw = Ji + Ja (3S) 



By use of the same series expansion as v/as previoiisly used, 
r 



^ 2Tr / -. m .^m 



-s + i 



sin m9dR 
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sin in 9 1 



2nir'^ m-s+2 



(r 



m-s+3 



- l) when + 



s in mS,, , ^ ^ 
>~**— - log r when 2n-'S + 2 = 0 

2mr^ 



(39) 



In the same manner 



J=> = - ~ / - rmrs-"jn: sin m e da 
2Tr 7 ^ Hin+s«i 

r 



sin m 9 



2in (m-.s-.2)r^^^^ 



(^•^o) 



From Gciuations (37). (3^), (39), and (Uo),.thG solu^ 
tion of equation (29) » '^ov the case in which 



is cons equent ly 

CD = si n ja^Q^ Fi s- g) 



V 



(when m-2?^s) 



cp = z^lJl^ J i + 1 log 
mr^ ^ ^ • 



2in 



(when in+2=s ) (Ul ) 



For the other cases the integration process follovrs 
exactly the same procedure. 
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Mach number when circulation is 0. 



NACA Fig. 2 



-4.3 



-4,2 



Glauert-Prandtl 
formula 

Kaman-Tsien 
formula 

One term approx- 
imation 

Two term. approx- 
imation 



0) 



?H 
0) 

o 
o 

<1> 
u 

01 




.2 .3 

Mach n'ombor, M 



Figure 2.- Variation of minimijun pressure coefficient 

with Maeb number when circulation is i/rU. 

4 



NACA 



Fig. 3 




.2 .3 

Mach nmber, M 



Figure 3.- Variation of miniraum pressure coefficient 

with Mach nxanber when circulation is ^nU. 
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